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We review different approaches to the graphical generation of the tadpole- free Feynman diagrams 
of the self-energy and the one-particle irreducible four-point function. These are needed for calcu- 
lating the critical exponents of the euclidean multicomponent scalar </> 4 -theory with renormalization 
techniques in d = 4 — e dimensions. 
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q ■ I. INTRODUCTION 

o : 

In 1982 Hagen Kleinert proposed a program for systematically constructing all Feynman diagrams of a field theory 
together with their proper weights by graphically solving a set of functional differential equations It relies on 
considering a Feynman diagram as a functional of its graphical elements, i.e., its lines and vertices. Functional 
derivatives with respect to these graphical elements are represented by removing lines or vertices of a Feynman 
diagram in all possible ways. With these graphical operations, the program proceeds in four steps. First, a nonlinear 
functional differential equation for the free energy is derived as a consequence of the field equations. Subsequently, 
this functional differential equation is converted to a recursion relation for the loop expansion coefficients of the free 
I ' energy. From its graphical solution, the connected vacuum diagrams are constructed. Finally, all diagrams of ri-point 
functions are obtained by removing lines or vertices from the connected vacuum diagrams. This program was recently 
£f) • used to systematically generate all Feynman diagrams of QED || and of 4 -theory in the disordered, symmetric phase 
ON I H and the ordered, broken-symmetry phase [@||. 

The present article reviews a modification of this program The aim is to construct directly the Feynman 
diagrams of n-point functions which are relevant for the renormalization of a field theory. To this end we consider the 
self-energy and the one-particle irreducible four-point function of the euclidean multicomponent scalar </> 4 -theory as 
functionals of the free correlation function. As such they obey a closed set of functional differential equations which 
can be turned into graphical recursion relations. These are solved order by order in the number of loops, producing 
all one-particle irreducible diagrams with their proper weights. A subsequent absorption of all tadpole corrections in 
the lines leads to modified graphical recursion relations for the tadpole-free one-particle irreducible diagrams which 
are needed for calculating the critical exponents with renormalization techniques in d = 4 — e dimensions. Finally, we 
elucidate how our procedure is related to the method of higher functional Legendre transformations which was also 
investigated in Ref. Q). 

^ ■ 

II. SCALAR -THEORY 

Consider a self-interacting scalar field </> with N components in d euclidean dimensions whose thermal fluctuations 
are controlled by the energy functional 



o 
o 



(2.1) 



In this short-hand notation, the spatial arguments and tensor indices of the field fa the bilocal kernel G 1 , and the 
quartic interaction V are indicated by simple number indices, i.e., 



l = {a; 1 ,a 1 }, / = / d d x 1 , fa = cj> ai (x x ) , G^ 1 = G Ql 1 Q2 (xi,x 2 ) , Vj. 

J 1 ™. « 



234 = Va t (xi,x 2 ,x 3 ,x 4 ) . (2.2) 



The kernel is a functional matrix G _1 , whi le t he interaction V is a functional tensor, both being symmetric in their 
respective indices. The energy functional ( |2.l| ) describes ci-dimensional euclidean </> 4 -theories generically. These are 
models for a family of universality classes of continuous phase transitions, such as the 0(A)-symmetric 4 -theory 
which serves to derive the critical phenomena in dilute polymer solutions (N — 0), Ising- and Heisenberg-like magnets 
(N = 1,3), and superffuids (N = 2). In all these cases, the energy functional (|2.ip is specified by the bilocal kernel 
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G-l a ,{x 1 ,xi)=6 au0t2 (-d 2 Xi +m 2 )5(x 1 -x 2 ), (2.3) 

and by the quartic interaction 

} 6{xi - x 2 )S(x! - x 3 )6(xi - X4) , (2.4) 
where the bare mass m 2 is proportional to the temperature distance from the cri tica l point, and g denotes the bare 



coupling constant. In this article we leave the kernel G 1 in the energy functional (2.1) completely general, except for 



the symmetry with respect to its indice s, a nd insert the physical value (2.3) only at the end. By doing so, we consider 
all statistical quantities derived from ( |2.l| ) as functionals of the free correlation function G which is the functional 
inverse of the kernel G _1 : 

G12 G 23 — 5\z ■ (2-5) 



2 

This allows us to introduce functional derivatives with respect to G whose basic rule reflects the symmetry of its 
indices: 

^ = Jtfl3*42+*uW. (2.6) 

dG 3 4 2 

Such functional derivatives are represented graphically by removing one line from a Feynman diagram in all possible 
ways ■ Thereby each line in a Feynman diagram represents a free correlation function 

1 2 = G 12 , (2.7) 

and each vertex represents an integral over the interaction 

V 123i • (2.8) 

1234 

Thus the differentiation (|2.6|) is illustrated graphically as 



2 = - I 1^34^2 + 1^43-^2 1, (2.9) 



S3 4 2 

where the elements of Feynman diagrams are extended by an open dot with two labeled line ends representing the 
delta function: 

1^2 = S 12 . (2.10) 



III. CONNECTED VACUUM DIAGRAMS 

By using natural units in which the Boltzmann constant kg times the temperature T equals unity, the partition 
function is determined as a functional integral over the Boltzmann weight e~ E ^\ i.e. 



Z = JV(f>e-* w , (3.1) 

and may be evaluated perturbatively as a power series in the interaction V. ^From this we obtain the negative free 
energy W = In Z as an expansion 

00 

W = ^W {1) , (3.2) 
1=1 

where the coeffi cien ts W^ 1 ' for each loo p ord er I > 2 may be displayed for as connected vacuum diagrams constructed 
from the lines (2.7) and the vertices (2.8). As has been elaborated in detail in Ref. |s|, the connected vaccum 
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diagrams contributing to follow together with their weights from a graphical recursion relation which can be 

written diagrammatically for I > 2 as 



W {1+1) = 



S 2 W^ 



1 
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6 I Si 2 S3 4 3 
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2 I Si 2 2 



This is iterated starting from the two-loop contribution 

wi2) = i CO 



l v — 

37 ~Si 



k=l 



- 2 2 



3 <5F^( fc+1 ) 

4 S3 4 



(3.3) 



(3.4) 



The right-hand side of ( [3.3] ) contains three different graphical operations. The first two are linear and involve one or 
two line amputations of the previous perturbative order. The third operation is nonlinear and mixes two different 
one-line amputations of lower orders.^] 

The connected vacuum diagrams resulting from the graphical recursion relation ( fr.3| ) toget her with their weights 
are shown up to I = 5 loops in Ref. ||. There we observed that the nonlinear operation in (3.2) does not lead to 
topologically new diagrams. It only corrects t he w eights generated from the first two linear operations. Continuing 
the solution of the graphical recursion relation (3.3) to higher loops is an arduous task. We have therefore automatized 
the procedure in Ref. || by computer algebra with the help of a unique matrix notation for Feynman diagrams. The 
corresponding MATHEMATIC A program and higher-order results up to I = 7 are available on the internet (7j . 



IV. ONE-PARTICLE IRREDUCIBLE DIAGRAMS 

Once the connected vacuum diagrams are known, the diagrams of the fully-interacting two-point function 



Gu = — I V> ><),,,_> < 



-E\, 



!>102 9394: e 



E 14>] _ Q^ 2 Q M — G 13 £?24 ~ G14G23 



and the connected four-point function 
G 1234 = \ J V< i> 

follow from removing one or two lines in all possible ways, respectively 

8W 



G12 = 2 1 G13G24 
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SG 



34 



°"1234 



5G\2 

2 / G35G46 -77: G13G24 — G14G 

0G56 



■■/■ 

J 56 



14^23 • 



(4.1) 
(4.2) 

(4.3) 
(4.4) 



Note that removing a line from the diagrams of G12 according to the first term of (4.4) leads to disconnected diagrams 



which are canceled by the second and the third term to yield the connected diagrams contributing to G\ 23i - Dropping 
all diagrams of Gn and £^234 which would fall into two pieces by removing a line, one obtains the one-particle 
irreducible diagrams of the self-energy 



S12 — G-^ 1 — G 12 



(4.5) 



with G 1 being the functional inverse of G 



1 Note that the first two operations in (3.3) were already used in Ref. |8| as heuristic algorithms to generate all topological 
different connected vacuum diagrams of the 4 -theory up to I = 8 loops with a computer program. Their corresponding weights 
were then determined by combinatorial means with a second computer program. 
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J G12 = S 
and the one-particle irreducible four-point function 



234 



5678 



— 1 s~i — 1 s~i — 1 s~i — 1 /"fC 

15 "26 "37 "48 "5678 • 



(4.6) 



(4.7) 



Along these lines all diagrams and their weights were constructed which are relevant for the five-loop renormalization 
of the <^> 4 -theory in d = 4 — e dimensions J8| — |l0|| ■ 

At this stage the question arises whether the one-particle irreducible diagrams can also be generated in a direct 
graphical way without any reference to the connected vacuum diagrams. This is indeed possible and will be elaborated 
in detail in Ref. ||. Here we restrict ourselves to present some of the results. To this end we extend the elements of 
Fcynman diagrams by a double straight line representing the fully-interacting two-point function 



2 = G 



12 



(4.8) 



and by a 2- and 4-vertex with a big open dot representing the self-energy and the one-particle irreducible four-point 
function 



2 3 
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1 4 



El2 , 
Tl234 



(4.9) 
(4.10) 



respectively. It turns out that the self-energy E follows from an integral equation which reads graphically 

Thereby the connected two-point function G is obtained from ( |4.5| ) according to the Dyson equation 

l — 2 = 1 2+1 == 2 • 



(4.11) 



(4.12) 



Solving the Eqs. (4.11) and (4.12) iteratively necessitates the knowledge of the one-particle irreducible four-point 
function T. Its diagrams could be determined from 

2 3 




(4.13) 

However, such a procedure would have one disadvantage: removing a line in the diagrams of the self- energ y E also 
leads to one-particle reducible diagrams which are later on cancelled by the third and the fourth term in ( 4.13 ) . As the 
number of undesired one-particle reducible diagrams occuring at an intermediate step of the calculation increases with 
the perturbative order, this procedure is quite inefficient in determining the diagrams of the one-particle irreducible 
four-point function T. By inserting (4.11) in ( 4.13] ) it turns out that we can derive another equation for T whose 
iterative solution only involves one-particle irreducible diagrams: 



1 4 

2 , 




(4.14) 



Note that Eqs. (4.11), ( 4.12| ) and (4.14) only generate one-particle irreducible diagrams for E and T, if all lower 
loop orders contain one-particle irreducible diagrams. By induction this establishes that E and T only consist of 
one-particle irreducible diagrams. 
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V. TADPOLE-FREE ONE-PARTICLE IRREDUCIBLE DIAGRAMS 



In order to reduce the number of diagrams, we aim at substituting the free correlation function G by a modified 
one G. If we would fix G according to 



G12 



G7 



(5.1) 



:S4 



it would contain all repetitive one-loop corrections. This method was established in Ref. [[Il] to get rid of all one- 
particle irreducible diagrams carrying tadpole corrections when calculating the f3- func tion of the vacuum energy 
density up to five loops. Here, however, we go one step further by demanding instead of (5.1) 



G"i2 — @12 + o / ^1234 G34 , 

1 J34 



(5.2) 



which amounts to absorbing all momentum- independent line corrections into the mass |T^ |. The tadpole corrections 
of the modified correlation function G arising from (5.2) were already treated perturbatively in Ref. It was shown 
that they lead to additional diagrams which cancel order by order all diagrams of one-particle irreducible n-point 
functions carrying any kind of tadpole correction. In Ref. || we elaborate that these inefficient cancellations can 
be circumvented as the remaining tadpole-free one-particle irreducible diagrams directly follow from a closed set of 
functional differential equations. To this end we consider the modified self-energy 



^12 



G\2 



(5.3) 



and the one-particle irreducible four-point function T as functionals of the modified correlation function G. Extending 
the elements of Feynman diagrams by a wiggly line representing the modified correlation function 



1 2 = Gx 2 

and by a 2-vertex with a wiggly dot representing the tadpole-free self-energy 

= £12 



our results read as follows. The tadpole-free self-energy £ obeys the integral equation 



(5.4) 



(5.5) 



(5.6) 



where the fully-interacting two-point function G c is obtained from (5.3) according to the modified Dyson-equation 

1^=2 = 1 2 + 1 — 0=2. (5.7) 

Then the one-particle irreducible four-point function T obeys 



2 3 

X 

1 4 



X 

1 < 

1 3 
2 , 



1 



>-X 4 

"-4 



7 Se~ 



1 4 

2 , 




1 2 

2 , 



(5.8) 



Iteratively solving Eqs. (5.6)— (5.8) in a graphical way leads to all tadpole-free diagrams of the self-energy and of the 
one-particle irreducible four-point function together with their weights. Up to five loops they are listed in Appendix 
A of Ref. Their respective quadratic and logarithmic divergencies in d = 4 — e dimensions contribute to the 

1 /e-poles of the renormalization constants of the field 0, the coupling constant g and the mass m 2 within the minimal 
subtraction scheme, so that they determine the critical exponents of scalar 4 -theory. 
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VI. ONE-PARTICLE IRREDUCIBLE DIAGRAMS WITHOUT LINE CORRECTIONS 



The number of diagrams can be even further reduced by substituting the free correlation function G by the fully- 
interacting two-point function G itself. Such a substitution was investigated in Ref. [j]] within the framework of 
higher functional Legendre transformations. By considering the self-energy E and the one-particle irreducible four- 
point function T as functional of the fully-interacting two-point function G, it can be shown that they obey pU|] 



2 3 
1 4 



1 

2 i = 
(5i 



53- 




Inscrting (3.1) in fl6.2|) it turns out that the one-particle irreducible four-point function also follows from 



(6.1) 
(6.2) 



2 3 
1 4 



1 1 

3 



5 ^X. 

6 5 ^ ^ 4 

7 5 6^=7 





The graphical solution of the functional differential equations ( J6 . l| ) and (6.2) or (6.3) leads to all one-particle irreducible 
diagrams which do not contain any line corrections. Once they are generated, we can recover the tadpole-free one- 
particle irreducible diagrams according to the following algorithm M. At first we subtract the one- loop correction 
from the self-energy in order to obtain the tadpole-free self-energy 



1 

2 i- 



(6.4) 



Iterating the modified Dyson equation (5.7), we then determine the tadpole- fre e dia grams of the fully-interacting two- 
point function G where the lines represent the modified correlation function (5.4 ). F inally we insert this expansion 
into the one-particle irreducible diagrams of E and T determined from (6.1) and (|6.2[) or fl6.3p. 



VII. SUMMARY 



We have reviewed different approaches to graphically generate tadpole-free one-particle irreducible diagrams to- 
gether with their weights which ar e ne e ded for calculating the critical exponents of <f> -theory. One approach is based 
on the graphi cal s olutio n of Eqs. (5.6)-( |5~8| ), the latter being more complex. The alternative approac h co nsists of the 
simpler Eqs. (|6.l|) and (|6.2[ ) at the expense of subsequent iterations of the modified Dyson equation (5/7) with ([Tl|). 
In order to decide which of those approaches is more efficient, it is necessary to perform further analytic studies or to 
automatize them by computer algebra. 
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